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ABSTRACT 

We present a study of correlations between spectral and timing parameters for a sample 
of black hole X-ray binary candidates. Data are taken from GX 339-4, H 1743-322, and 
XTE J1650-500, as the Rossi X-ray Timing Explorer (RXTE) observed complete outbursts 
of these sources. In our study we investigate outbursts that happened before the end of 2009 
to make use of the high-energy coverage of the HEXTE detector and select observations that 
show a certain type of quasi-periodic oscillations (type-C QPOs). The spectral parameters are 
derived using the empirical convolution model simpl to model the Comptonized component 
of the emission together with a disc blackbody for the emission of the accretion disc. Ad- 
ditional spectral features, namely a reflection component, a high-energy cut-off, and excess 
emission at 6.4 keV, are taken into account. Our investigations confirm the known positive 
correlation between photon index and centroid frequency of the QPOs and reveal an anti- 
correlation between the fraction of up-scattered photons and the QPO frequency. We show 
that both correlations behave as expected in the "sombrero" geometry. Furthermore, we find 
that during outburst decay the correlation between photon index and QPO frequency follow a 
general track, independent of individual outbursts. 

Key words: X-rays: binaries - X-rays: individual: GX 339-4, H 1743-322, XTE J1650-500 
- binaries: close - black hole physics 



1 INTRODUCTION 

Most known black hole X-ray binaries (BHT) are transient. 
They are only observable during outbursts, as they are too 
faint to be detectabl e with present X- ray instruments during 
quiescence (see e.g. iGarcia etak I Il998h . While a BHT is in 
outburst it evolves through different states, which show char- 
acteristic timing and spectral properties. The states can be 



identified with the help 


of the hardness intensity diagram 


(HID; Homanetal. 200 ll: 


Belloni etal, 20051; Homan & Bellonil 


20051; Gierlinski & Newton 


20061; McClintock & Remillard 20061; 


Fender et alj 2009; Belloni et alj|201 ll). the hardness-rms diasram 


(HRD; BeUoni et al. 200^. and the rms-intensitv diagram (RID; 



the low hard state (LHS) and there is a transition to the high soft 
state (HSS) in between. Three different types of low frequency 
quasi-periodic oscillations (LFQPOs) can be distinguished in BHTs 
jCasella et al.l2005l ; IWiinands et al.ll99^) . In the LHS and the hard 
interm ediate state (HIMS; following the nomenclature of iBelloriil 
I2OIOI) a specific timing fe ature named type-C QPOs can be ob- 
served l lBelloni et al J201 ih . 



* E-mail: holger.stiele@brera.inaf.it 



Existence of a correlation b etween the QPO frequency a nd the 
photon index was first shown bv lDi Matteo & Psaltij l ll999l) . They 
focused their discussion on the implications of the limited QPO 
frequency b and on the limited change in the irm er disc radius during 
transitions. S haposhnikov & TitarchukI ( l2009t) studied correlations 
between the centroid frequency of QPOs and spectral parameters 
for a sample of eight BHTs observed with the Rossi X-ray Timing 
Explorer (RXTE). To obtain the spectral parameters RXTE/PCA 
(Proportional Counter Array) spe ctra in the 3.0 - 50.0 keV range 
were fitted with the bmc model jTitarchuket"^! 19971) , which is 
hardwired to a Planck function. 



In this paper we investigate correlations between spectral pa- 
rameters and timing properties for a sample of seven outbursts from 
three different sources. The spectra are fitted with a disk blackbody 
convolved with the simpl model (Steiner et al. 2009), which is an 
empirical model of Comptonization. To model the broad spectral 
features due to reflection accurately, and to determine a high-energy 
cut off - if present -, we included RXTE/HEXTE data in our study. 
We used only observations that showed type-C QPOs. A discussion 
of the correlation between photon index and centroid frequency for 
type-B QPOs observed in the 2010 outburst of GX 339-4 can be 
found in lStiele et all ( I2OIII) . 
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Figure 1. Long-term PCA light curves of GX 339-4 (upper panel), 
XTE J1650-500 (middle panel), and H 1743-322 (lower panel). Each point 
corresponds to an entire observation. The observations used in our study are 
marked with different symbols. 



2 OBSERVATIONS AND DATA ANALYSIS 

In this paper we analysed archival RXTE observations of several 
outbursts of the black hole X-ray binary candidates GX 339-4, 
H 1743-322, and XTE J1650-500, namely: 

• the 2002/03, 2004, and 2007 outbursts of GX 339-4 

• the 2003, 2008, and 2009 outbursts of H 1743-322 

• and the 2001 outburst of XTE J 1650-500, which is the only 
observed outburst of this source. 

We selected these sources as RXTE observed complete out- 
bursts of them, i. e. neither the initial rise nor the decay to quies- 
cence is missing; a substantial number of observations with type- 
C QPOs have been detected; and the HIDs were nicely q-shaped. 
In addition, RXTE observed several outbursts of GX 339-4 and 
H 1743-322, including outbursts at different luminosities and one 
outburst of H 1743-322 that did not go all the way to the soft state 



jMotta etal.l201(]| . l201lh . Long-term PCA light curves are given in 
Fig. [T] Outbursts that happened either partially or totally after the 
end of 2009 have been excluded from our study to avoid the intro- 
duction of systematic uncertainties in the spectral parameters due to 
the strong residuals in the HEXTE spectra that are related to the dif- 
ficulties in determining the background contribution s ince techni- 
cal pr oblems occurred in the HEXTE detector (see also lStiele et al.l 
l2012h . We investigated the timing and spectral properties, as de- 
scribed below. 



2.1 Timing analysis 

The values of the centroid frequencies and the classification of 
th e QPOs for a l l three outbursts of GX 339-4 have been published 
in iMotta et all (|20 1 ll), and for th e 2008 and 2009 outbursts of 
H 1743-322 in iMotta et al] 1 I2OIOI) . The data of the 2001 outburst 
of XTE J 1650-500 were analysed using the methods and clas- 
sification criteria described in these papers. We used data from 
the Proportional Counter Array (PCA) to compute power density 
spectra (PDS) for each obs ervation following the procedure out- 
lined in lBelloni et al. liooi). PDS production has been limited to 
the PCA channel band - 35 (2 - 15 keV) and used 16 sec- 
ond long stretches of Event mode data. For observations which 
showed a quasi-periodic oscillatio n (QPO) we subtra cted the con- 
tribution due to Poissonian noise llZhangetal.lll995l) . normalised 
the PDS accor ding to Leahy et al. ( 198^ and converted it to square 
fractional rms l lBelloni & IIasinger 19900 . We determined the cen- 
troid QPO frequency, by fitting the noise c omponents as well a s 
the QPO feature with Lorentzians, following iBelloni et al.l ( l2002h . 
PDS fitting was ca rried out within the standard xspec fitting pack- 
age jArnaud|[l996l) by using a one-to-one energy-frequency con- 
version and a unit response. The observations of X TE J 1650-500 
during outburst decay have been already analysed in lKalemci et al.l 
( I2OO3I) . They used a 2 - 30 keV energy range with 128 second long 



stretchefl For the 2003 outburst of H 1743-322 we t ook th e QPO 
frequencies listed in Table 2A of iMcClintock et al] l l2009l) . They 
used the full bandwidth of the PCA instrument (2 - 40 keV) and 
searched for QPOs in the 4 mHz to 4 kHz range w ith a sliding 
frequency window technique jRemillard et alj|2002al lbl). In the fol- 
lowing, we focus on observations which QPOs we have classified 
as type-C. 



2.2 Spectral analysis 

For all observations with type-C QPOs we used the PCA Standard 
2 mode (STD2), which covers the 2-60 keV energy range with 129 
channels, and the HEXTE Standard mode, which covers the 15- 
250 keV energy range with 129 channels, for the spectral analysis. 
The standard RXTE software within heasoft V. 6.9 was used to 
extract background and dead-time corrected energy spectra for each 
observation. Solely Proportional Counter Unit 2 from the PCA was 
used since only this unit was on during all the observations. To 
account for residual uncertainties in the instrument calibration a 
systematic error of 0.6 per cent was added to the PCA spectral. For 
observations taken between 2001 and 2004 we used HEXTE data 



' Using 16 second long or 128 second long stretches does not lead to any 
change in the derived QPO centroid frequency, as type-C QPOs aie strong, 
nanow features in the PDS. 

^ A detailed discussion on PCA calibration issues can be found at: 
|http://www. universe.nasa.gov/xrays/programs/rxte/pca/doc/nnf/pcarmf- 1 1 .7/ 1 



© 2012 RAS, MNRAS 000,[T]|8] 



60113-01-03-00 



Relations between timing and spectral properties 

60113-01-13-02 




10 20 



50 100 



10 20 



50 100 



Energy (keV) 
60118-02-11-00 



Energy (keV) 
60118-02-15-00 




H \ — I — 



Figure 2. 

the decay 



5 10 20 50 100 

Energy (keV) 
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from cluster A, while for observations taken in 2007, 2008, and 
2009 cluster B data have been used. 

Co mbined PCA+HEXTE sp ectra were fitted within isis 
V. 1.6.1 l lHouck & Denicolal 120001) in the 4-40 keV and 22 - 
200 keV range. We uniformly fitted the spectra with a partially 
Comptonized multi-color disc blackbody model, including fore- 
ground absorption. Four vF,, spectra of XTE J 1650-500 taken form 
the beginning and end of the rise and decay branch, respectively, 
are shown as exemples in Fig. [2] isis "unfolded spectra" are in- 
dependent of the assumed spectral model (i.e., the unfolding is 
done solely with t he response matrix and effective area files; see 
iNowak et al.ll200 j for details). The plotted residuals, however, are 
those obtained from a proper forward-fol ded fit. The disc emi ssion 
was approximated by the diskbb mo del jMitsuda etaLlll984h and 
the simpl model jSteiner et al.ll2009[) was used for Compton scat- 
tering. The latter one being an empirical convolution model that 
converts a given fraction of the incident spectrum into a power law 
shape with a photon index F. The amount of the up-scattered frac- 
tion of the incident radiation is stored in a parameter called scat- 
tered fraction. We allowed for a reflection component and mod- 



eled it with reflect jMagdziarz & Zdziarskilll995h . The reflec- 
tion component aff'ected solely the up-scattered photons (see also 
IGou et al.ll201ll ; iTamura et alj|2012h . As simpl and reflect are 
convolution models it is necessary to calculate the model well- 
outside the normal bounds of the PCA and HEXTE energy range. 
Therefore, we used an extended energy range from 0.1 keV to 1 
MeV for our fits. The distribution of reflection values for all three 
sources are shown in Fig.[3l The mean error is about 0. 15 (GX 339- 
4, XTE J1650-500) to 0.2 (H 1743-322), although individual obser- 
vations can show much smaller or bigger errors. If needed a high- 
energy cut-off (highecut) was included. This was mostly the case 
for observations observed during the rise of an outburst. We added 
a Gaussian to account for excess emission at 6.4 keV. The centroid 
was allowed to vary between 6.4 and 6.7 keV and the line width 
was constrained between and 1 keV to prevent artificial broaden- 
ing due to the response of the PCA detector at 6 .4 keV. For the fore - 
ground absorption we used the TBabs model l lWilms et alJbOOOl) . 
with fixed foreground absorption (A'h)- Values us ed for the indi- 
vidual sources are: 5.0x10^' cm - for GX 339-4 |Zdziarski et"al] 
[2004i) , 1.6x10-^ cm-- for H 1743-322 jCapitanio et al.ll2009l) , and 
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Figure 3. Distribution of the reflection values for GX 339-4 (upper panel), 
XTE J 1650-500 (middle panel), and H 1743-322 (lower panel) with a bin- 
ning of 0. 1. The mean error is about 0. 15 (GX 339-4, XTE J 1650-500) to 0.2 
(H 1743-322). Open bars denote observations belonging to the rise branch 
of an outburst, while filled bars mark observations of the decay branch. Dif- 
ferent outbursts are marked with dift'erent colours (grey shades), using the 
same colouring scheme as in Figs.[T]and|4] 



6.7x10-' forXTE J1650-500 fromsick et alj2004l) . For most 
observations the disc blackbody temperature lies between 0.5 and 
1 .2 keV and we obtain an inner disc radius of a few tens kilometers. 
The reduced value, x^^i = I^Aof, where N^^f is the number of 
degrees of freedom (dof), is less than or around one for most of the 
observations. For a small fraction (« 2%) of spectra the value of 
xl^^ exceeds 1.5. However, it never reaches a rejection limit of two. 



3 CORRELATIONS 

In this section we present correlations between QPO centroid fre- 
quency and different spectral parameters. Figure |4] shows for the 
three BHTs the correlations between QPO centroid frequency and 
photon index or scattered fraction, respectively. Individual out- 
bursts of the same source are marked with different symbols. The 
broad distribution of reflection fractions (see Fig. [3} and the rather 
large systematic errors hamper the detection of any clear trends 



in the correlation between reflection fraction and QPO frequency. 
We do not show correlations between QPO frequencies and disc 
parameters, as the obtained disc parameters should be taken with 
care. The working range of PCA (3 - 40 keV) covers only the 
high energy part of the disc component, above the Wien peak. 
The missing coverage of lower energies together with the pres- 
ence of a strong Comptonized component increases the uncertain- 
ties in the derived disc parameters, especially for observations ob- 
tained at the very begin and end of an outburst. In addition, it is 
known that the spectral parameters derived fro m the diskbb model 
should not be interpreted liter ally (see e. g. iMerloni et alj|200^ : 
iRemillard & McClintockll200d) . Nevertheless, we want to mention 
that our investigations imply that QPOs at the same frequency oc- 
cur at a lower inner disc radius in the decay branch compared to the 
rise branch. 



3.1 The scattered fraction - QPO centroid frequency relation 
(SF-QPO relation) 

All investigated outbursts show a negative correlation between the 
fraction of up-scattered photons and the QPO centroid frequency: a 
high scattered fraction is observed at low QPO frequencies, and the 
scattered fraction decreases with increasing frequency. In the case 
of XTE J 1650-500 the correlation is most obvious and we recog- 
nise two branches which belong to the rise and decay branch of 
the outburst, where the latter one lies above the branch of the rise. 
For the other two sources the distinction between rise and decay 
is less clear. For H 1743-322 it is impossible to separate different 
branches. For a few observations we obtained a scattered fraction 
close to unity. In these cases we used the Icr lower error as a lower 
limit (indicated in Fig.|4]by a small arrow). 



3.2 



The photon index 
(F-QPO relation) 



- QPO centroid frequency relation 



The photon index rises with increasing QPO centroid frequency. In 
the case of XTE 11650-500 it is a shallow and rather linear rise. The 
two branches belonging to the rise and decay of the outburst are in- 
distinguishable. We recognise distinct rise and decay branches for 
the 2002 and 2007 outburst of GX 339-4. The branch belonging 
to the rise of the 2004 outburst seems to be aligned with the de- 
cay branch of the 2002 outburst. This might be related to the lower 
luminosity at which the 2004 outburst has been observed. The pho- 
ton indices are well constrained. The F-QPO relation in GX 339-4 
can be described as a loop-like structure compared to the narrower, 
more band-like appearance in the case of H 1743-322. 



3.3 F-QPO relation during outburst decay 

In Fig.|5] (left column) we summarised in one diagram the F-QPO 
relation during decay for all outbursts. This figure reveals that the 
F-QPO relation follows a general track during outburst decay. The 
relation can be approximated by F = Fq -H a * cfgpo, with Fq = 1.73 
and a = 0.08. The standard deviation in F is 0. 10. We would like to 
point out that during outburst decay also the distribution of source 
luminosities is narrower. 



3.4 SF-QPO relation during outburst decay 

In Fig.[5] (right column) we summarised in one diagram the SF- 
QPO relation during decay for all outbursts. Unlike the F-QPO rela- 
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Figure 4. F-QPO (left column) and SF-QPO (right column) correlations for GX 339-4 (upper row), XTE 11650-500 (middle row), and H 1743-322 (lower 
row). Open symbols denote observations belonging to the rise branch of an outburst, while filled symbols mark observations of the decay branch. If a source 
has been detected during more than one outburst, the individual outbursts are marked with dift'erent symbols. Observations for which we obtained only a lower 
limit of the scattered fraction are marked by arrows. 



tion, the SF-QPO relation shows two tracks. The upper one consists 
of all observations of XTE J 1650-500 and of most observations of 
H 1743-322. The lower track contains the observations of GX 339- 
4. Only, the 2003 outburst of H 1743-322 starts at ~7 Hz on the cor- 
relation of GX 339-4, leaves it and reaches the upper track spanned 
by the other observations of H 1743-322. 



4 DISCUSSION 

We studied for the BHTs GX 339-4, H 1743-322, and XTE J1650- 
500 the correlation between the QPO centroid frequency and two 
spectral parameters: the photon index and the fraction of up- 
scattered photons. The correlations presented in this work agree 
qualitatively with the ones that can be deri v ed fro m the values given 
in Table 2 of IShaposhnikov & TitarchukI ( |2009|) . This means that 
the anti-correlation in the SF-QPO relation and the positive correla- 
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Figure 5. The F-QPO (left column) and SF-QPO (right column) relation during decay for all seven outbursts. In the F-QPO diagram all points seem to follow 
the same track, without recognisable differences between individual sources. In the SF-QPO diagram there are two different tracks. The points obtained from 
observations of H 1743-322 and XTE J1650-500 are located on one track, while those of GX 339-4 span their own track at lower scattered fractions. 



tion in the F-QPO relation are model independent. We note that the 
photon indices in the present work are slightly higher than those de- 
rived with the bmc model. This has to be exp ected, as we inc l uded a 
reflection component in our spectral model. IVignarca et Zl | |2003|) 
investigated the correlation between QPO frequency and photon in- 
dex for an additional five sources, partially taking QPO frequencies 
and spectral parameters from the literature. A single correlation be- 
tw een spectral inde x and QPO centroid frequencies was also found 
bv iKalemcil ( l2002h during the outburst decays of six sources in 
eight outbursts. 

X-ray spectra of black hole X-ray binaries consist of two main 
components: The soft component is believed to originate in the geo- 
metrically thin and opticall y thick Shakura-Sunyaev accretion disc 
jShakura & Sunvaevlll973l) . One of the most plausible processes 
of formation of the hard spectral co mponent is Comptonization 
of soft disc photons on hot electrons jSunvaev & Truempeilll97^ ; 
ISunvaev & TitarchukllT980l) . The Comptonization site is often re- 
ferred to as a corona. Although it is generally accepted that the 
Comptonizing corona has to be located in the close vicinity of the 
black hole, there are currently different ideas brought up on the 
detailed geometry of the region. The overall behaviour of the SF- 
QPO and F-QPO relations can be explained qualitatively withi n 
the "sombrero" geometry jPoutanen et"al]|l997l ; [Gilfanovll2010h . 
In this configuration, a quasi-spherical corona surrounds the black 
hole and the accretion disc extends a short distance into the corona. 
In the LHS the disc is truncated at a large radius and the remaining 



space is filled by the hot, optically thin corona llDone et al 



Gierliiiski et al. 



20081 ; iD'Angelo et"ai]|2008l : iTomsick et a~ 



2007 



200S 



Hiemstra et alj|2009h . While during the HSS, the disc reaches the 



inner most stable orbit leaving (at most) only little space for the 
corona. A fraction of the Comptonized photons emitted from the 
corona irradiates the accret ion disc. Some o f them are reflected 
due to Compton Scattering l lBaskoetal.lll974h . We note that there 
are claim s that the disc reach es the inner most sta ble orbit even in 
the LHS jMiller et alj[2006al lbl: lRvkoflFet alj|2007h . There are also 
models which predict a hot inner and a cool outer disc se parated by 
a gap filled with an advection-dominated accretion flow jLiu et al] 
l2007h . 

Coming back to the sombrero geometry, which matches our 
findings, the QPOs correspond to oscillations in a tra nsition layer 
between the disc and the hotter Comptonizing region l lNobili et al.l 



I2OOOI : llngram et al.l |2009| ; llngram & Pond [201lh . The QPO fre- 
quency is inversely related to the truncation radius of the disc , 
as it is in most models (see e . g. iTitarchuk & OsherovichI I2OO0I : 
IChakrabarti & Manickan]l2000l) . While the system evolves from 
the LHS to the HIMS the truncation radius moves inward and the 
QPO frequency increases. In the sombrero configuration the solid 
angle of the disc seen by the hot electrons correlate with the frac- 
tion of disc emission reaching the corona. As the disc expands to- 
wards the black hole, the reflection scaling factor increases, the 
frac tion of up-scattered photons decrease s and the spectrum steep- 
ens l lZdziarski et alJl999l : lGirfanovl2010l) . Thus the observed (anti- 
)correlations between QPO frequency and spectral parameters can 
be explained within the sombrero geometry. 



5 CONCLUSION 

We studied correlations between spectral and timing parameters for 
observations with type-C QPOs. The sample comprised observa- 
tions of GX 339-4, H 1743-322, and XTE J 1650-500. Our investi- 
gations confirmed the known positive correlation between photon 
index and centroid frequency of the QPOs and revealed an anti- 
correlation between the fraction of up-scattered photons and the 
QPO frequency. We showed that both correlations behaved as ex- 
pected in the "sombrero" geometry, which also predicted the ob- 
ser ved correlation b etween photon index and reflection scaling fac- 
tor ( lGilfanovll2010l) . Furthermore, we showed that during outburst 
decay the correlation between photon index and QPO frequency 
followed a general track, independent of individual outbursts. 
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